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We investigate the physical consequences of orbital current patterns (OCP) in doped two-leg Cu-0 
Hubbard ladders. The internal symmetry of the pattern, in the case of the ladder structure, differs 
slightly from that suggested so far for cuprates. We focus on this OCP and look for measurable 
signatures of its existence. We compute the magnetic field produced by the OCP at each lattice 
site, and estimate its value in view of a possible experimental detection. Using a renormalization 
group (RG) analysis, we determine the changes that are caused by the SU(2) spin-rotational sym- 
metry breaking which occurs when the OCP is present in the ground state phase diagram. The 
most signifcant one is an in-plane SDW gap opening in an otherwise critical phase, at intermediate 
dopings. We estimate the value of this gap, give an analytic expression for the correlation functions 
and examine some of the magnetic properties of this new phase which can be revealed in measure- 
ments. We compute the conductance in the presence of a single impurity, using an RG analysis. 
A discussion of the various sources of SU(2) symmetry breaking underscores the specificity of the 
OCP induced effects. 



I. INTRODUCTION 

More than two decades after the experimental discov- 
ery of a pseudogap in the phase diagram of high tem- 
perature superconductors (HTSC)», a microscopic expla- 
nation of its origin still remains elusive. Many scenar- 
ios have been proposed 2 - - — but so far no consensus has 
emerged on any particular one. Among the contenders 
are proposals that the pseudogap might be linked to cur- 
rent formation. It was suggested that staggered current 
patterns (DDW phase}^ would lead to a pseudogap due 
to the doubling of the unit cell. An alternative uniform 
current patter n 10 ' 11 was proposed by Varma. In this lat- 
ter case the three band nature of the system is crucial, 
since currents flow along closed loops formed by Cu-0 
bonds in such a way that translational symmetry is pre- 
served, while time-reversal symmetry is broken. On the 
experimental front specific signatures compatible with 
uniform currents where observed^ 2 - - — even if the issue 
is still controversial 16 ' 17 . On the theoretical side signifi- 
cant efforts were made to probe the possible existence of 
such a phase, either using numerical approaches^— or 
starting from Mott insulating side2i. 

As a way of investigating the physics of these states in 
a more controlled way and also in view of possible con- 
nections with experimental systems 2 ^—, phases exhibit- 
ing OCP were also considered in the context of two-leg 
ladder s 27 ' 28 . Indeed for this essentially one-dimensional 
(ID) case, various analytical tools allow one to explore 
the consequences of such an unusual phase. In previ- 
ous wor k 29 ' 30 we established that in the weak coupling 
limit of Cu-O Hubbard ladders there exists a range of 
doping with a quasi-long range ordering of orbital cur- 
rents, such that the pattern is symmetric with respect to 
the exchange of the two legs (o-OCP). In these studies, 
however, possible changes in the symmetry of the system 
caused by the formation of the orbital currents was not 
taken into account. Our analysis of the two-leg ladder 



problem was based on the assumption that the spin de- 
grees of freedom were totally decoupled from the motions 
of the carriers in real space. Although this is a very good 
starting point when the order parameter is small, it is 
interesting to investigate the consequences of such a cou- 
pling. In particular in the presence of a static pattern 
of currents, it is natural to expect that the spin of the 
electrons will couple to the generated moment, and lead 
to interesting spin-orbit effects. 

The aim of the present paper is to explore and to 
discuss physical effects which are occurring in response 
to the o-OCP. Anticipating the results presented below, 
the main consequence of the o-OCP is the appearance 
of magnetic moments which break SU(2) spin-rotational 
symmetry. While measuring the primary (first order) re- 
sponse of the system to the o-OCP is expected to be 
experimentally challenging, these collateral effects con- 
stitute clear fingerprints of the existence of such a phase. 
We also check the stability of the o-OCP itself with re- 
spect to the second order perturbation which couples the 
order parameter with the spin of the electrons. We ar- 
gue that SU(2) symmetry is fairly well protected (in the 
absence of crystalline anisotropies) and that its breaking 
is a telltale of the existence of the o-OCP. 

The paper is organized as follows: At the beginning 
of Sec. |H] we introduce the Hubbard model for two-leg 
Cu-O ladders and we review the main results that were 
obtained for this model in the weak coupling renormaliza- 
tion group (RG) approach. Next we derive the additional 
terms that arise in the presence of the o-OCP. The mag- 
netic moments that they induce in the ladder are calcu- 
lated and the first order effect that they generate are also 
discussed there. New terms in the hamiltonian will be de- 
rived. Sec. IIIII is devoted to an RG analysis of the gen- 
eralized ladder model obtained in the previous section. 
Experimental consequences entailed by these results are 
described in Sec. IIV1 Lastly, in the Appendix, we dis- 
cuss other possible sources of SU(2) symmetry breaking 
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in two-leg Cu-O ladders (such as Rashba effect). 

II. THE MODEL 

Let us begin with a brief review of two-leg Cu-O ladder 
physics. We consider a two-leg Cu-O ladder with local 
interactions. The Hamiltonian of this system contains 
two parts: the kinetic energy of electrons moving on the 
lattice Ht and electron interactions Hint 

H = H T + H mt (1) 

In the SU(2) invariant case, the physics of this hamilto- 
nian is known (assuming that Hi nt is treated as a per- 
turbation of the kinetic term) 29 ' 30 . In the first part of 



this section, we recall its main characteristics, we show 
the condition required for the o-OCP to exist, and we 
also introduce relevant notations. In the second part, we 
explain how the o-OCP breaks spin-rotational symmetry 
and we derive new terms in the hamiltonian that arise be- 
cause of the broken symmetry. In principle, these terms 
affect both Ht and Hint ■ This secondary effect is caused 
by the presence of the o-OCP. 



A. The Cu-O Hubbard ladder 

Including two Cu and five O atoms in the unit cell, the 
tight-binding kinetic energy part Ht reads 
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where 



is the creation operator of holes with 



spin a on a copper (oxygen) site (J is a site along chain 
and m labels the atoms in each cell); n ( j n n - a = a} m - a ■ a m j a 
We use hole notation such that t, t±,t pp are all positive. 
e = to ~ e Cu is the difference between the oxygen and 
copper on-site energies. 

The model is reduced to two low lying bands crossing 
the Fermi energy. They are denoted o (symmetric under 
the exchange of the two legs) and 7r (antisymmetric under 
the exchange of the two legs). With a (= 0,tt) and a 
denoting the band and the spin index respectively, the 
Hamiltonian reads 



H T = y^e a (fc)n fcc 



where 



&mko 



(3) 



(4) 



and e a (k) and A mQ are the eigenvalues and components of 
the eigenvectors of the Hamiltonian matrix (see Ref. [29h . 

In the low energy limit one may linearize the dispersion 
relation in the vicinity of the Fermi energy: 



Ht = r 1 V F a ^ arqa a a 

\q\<Q ra ° 



rqa 



(5) 



and it is easy to bosonize this free fermion theor y 31 ] 32 . 
Two charge (c) and spin (s) boson phase helds <j)(x) are 



introduced for each fermion species (x is the spatial co- 
ordinate along the ladder). 

We also introduce the phase fields 9; their spatial 
derivative 11(2;) = 7rV#(x) is canonically conjugated to 
4>(x). Now the Hamiltonian may be rewritten using the 
above phase fields. The kinetic part and those pieces 
in Hi n t which can be expressed as density-density terms 
give rise to the following quadratic form 



^ = E/|[(^)( 7 rn A ) 2 + (^ 



){dAxf] (6) 



in the diagonal basis Bq. X = 1,2,3,4 labels the eigen- 
modes (1,2 are spin modes and 3, 4 are charge modes^ 9 -) 
For the non-interacting system one has K\ = 1 for 
all modes; the diagonal density basis is then the bond- 
ing/antibonding one B 07T (the momentum kj_ associated 
with the rungs is either or tt) . The other basis which is 
commonly used in the literature is the total/transverse 
one, It is related to B ow by: 



W(-) 



V2 



(7) 



where v =c or s, depending on which particular density 
one considers. 
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The interaction part, in fermionic language, is given 

by 

+ E Uon m jt n mji 

j m£Cu m£0 

+ E y^.Vcu-Onmjannjv') (8) 

Eq. © gives rise to two types of terms: the first ones are 
of the forward scattering type and they can be cast in 
a quadratic form (~ V</>(x)) 2 in bosonization language. 
These terms are then incorporated in the above men- 
tioned matrix K of Luttinger liquid (LL) parameters, and 
hence they are treated exactly in this procedure. Since 
the form of K depends on the basis in which the densities 
are expressed, the Hamiltonian will take the simple form 
Eq. (J5]) in the eigenbasis of the matrix K. The remain- 
ing interaction terms yield non-linear cosines and these 
are the ones for which the RG procedure is required. In 
Refs. [29l|30l we showed that the eigenbases for the spin 
and for the charge modes rotate during the RG flow. Two 

fixed points were found, namely and B 01T . Interband 

physics dominates in the former case (mixing of the o and 
7r bands) and intraband physics, in the latter case. Close 

to half filling (low doping regime) , is the fixed point 

basis for the spin and the charge modes. Both spin modes 
are gapped, and so is one of the charge modes. Using the 
notation of Balents and Fisher- 3 the ladder is in the C1S0 
state (CnSm denotes n (m) gapless charge (C) (spin S) 

modes). For an intermediate range of dopings, B^ is the 

fixed point basis for the spin variables and i? 07r the fixed 
point basis for the charge variables. This is the C2S2 
regime when all modes are gapless. The decoupling of 
the spin and of the charge eigenbasis is responsible for 
this quantum critical state. For higher dopings, B ow is 
the fixed point basis for the spin and the charge modes. 
The ladder is in the C2S1 phase where the o spin mode 
is gapped. The corresponding phase diagram, which was 
established in Refs. I291I301 is shown in FigHJ 

A noteworthy result that came out of our analysis was 
that the symmetry of the OCP is that of the o band, 
for the Cu-O ladder. What is more, there exists a range 
of values of t pp , above some critical hopping where this 
phase is dominant. In the following, we present a detailed 
discussion of experimentally accessible effects caused by 
the presence of this unusual phase. 

B. Primary effect induced by the o-OCP: magnetic 
moments 

When the o-OCP is present, its current pattern I(x) 
generates a local magnetic field B(x) at each site x of 
the lattice. This is the quantity that we are interested 
in. The value of the field can be obtained from the Biot- 
Savart law 

ot \ fo / rM (Axr Mo t I' dlxr m\ 
B(x) = _^j(x)_ = _£j B £_ (9) 




FIG. 1: The phase diagram of two-leg Cu-O Hubbard ladders 
versus doping for Ucu > 0. Zero doping corresponds to the 
half filled case. Umklapp terms which open a gap in the charge 
symmetric mode are not included here^. SCd, CDW and 
SDW denote respectively d-wave superconductivity, charge- 
density wave and spin-density wave states. OCP indicates 
that an orbital current pattern is present on top of a density 
wave instability 



r is the vector joining the center of an infinitesimal ele- 
ment of current Idl to the lattice site x. In the second 
equality, the prime on the integral means that we are 
integrating over an elementary Cu-O loop (triangle) - 
denoted by n - on the ladder; we sum over all triangles 
and we use the fact that current is conserved for each 
loop n of the lattice. 

To perform the summation in Eq. ([9]), we first single 
out those triangles that belong to the particular unit cell 
that contains the lattice site x; we immediately notice 
that we get a zero net contribution whenever 

• the current flows in a direction that passes through 
the lattice site x (the cross product dl x r gives 
zero) 

• the high symmetry of the OCP causes a cancel- 
lation of the fields due to currents connected by 
mirror symmetry (for that reason, the field on the 
on-rung oxygen is equal to zero) 

These properties can easily be accounted for if one 
rewrites the magnetic field using a multipole expansion 
of the vector potential 

A W = £E 7 «E^T / (rTP m (cos6)dl (10) 

There are no magnetic monopoles, dipole contributions 
vanish because of the cancellations described above and 
so do quadrupole contributions inside each cell. Symme- 
try allows us to conclude that octupoles are the lowest 
order non- vanishing terms. The above formula can be 
significantly simplified thanks to the fact that all atoms 
lie in the same plane from which we have cos# = 1; we 
can also assume that the approximate distance between 
the moments r' is equal to unit cell size a. 
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Thus, we notice that we may get an accurate estimate 
of the contribution of all the other cells, if we assimilate 
them to octupole moments. Basically this approxima- 
tion amounts to replacing the primed integrals within 
each triangle by a number (the magnetic moment at the 
center of each elementary triangle). The current creates 
a dipole moment perpendicular to the plane of the trian- 
gle (and proportional to the current flowing around it). 
Four neighboring current triangles, centered at a single, 
n-th Cu site, define an elementary octupole Ci n ) (with a 
shape of a square). Then the magnetic field is approxi- 
mated as a sum of octupoles along the ladder: 

= A^o rCW , x^ Xjn) - ((-n) COO - C(-") u 
4ir l a 5 ^ [ r(x n ) 5 (r{x n ) 2 + a 2 ) 5 / 2 h 

The first term comes from octupoles on opposite legs, 
the sum runs along the rungs of the ladder, with the first 
piece arising from same leg contributions and the second 
piece from opposite leg contributions. A crucial point 
is how to estimate the value of £(n) ~ /„. This work 
pertains to physical systems quite similar to cuprates. 
Experimental magnitudes of the local magnetic (dipolar) 
moments, which are tilted 45° out of the Cu-O plane, 
are of order O.Ihet^. Thus we assume that our maximal 
octupolar momentum (on top of the density wave (DW)) 

is equal to C{I ma x) — 4^0.1/isa 2 and we will use that 
value in our estimates. The case of on-leg oxygens, which 
lie at the boundary of each cell, is special because then 
the quadrupolar contribution does not vanish, and it will 
be discussed below. 

One has to remember that the periodicity along the 
ladder plays an important role in the computation of such 
a sum; specifically, if the o-OCP were uniform along the 
ladder, all the different contributions would cancel out, 
because then £(n) = £(— n), and the total magnetic field 
would be zero on each atom. In our case we have an 
o-OCP on top of a DW with a real space periodicity 
(2fcj? ) _1 , hence a magnetic field with the same period- 
icity appears. 

As a result one finds that, depending on the nature of 
the lattice site x, the magnetic field contains either only 
octupolar or both quadrupolar and octupolar contribu- 
tions. 

• for each on-leg oxygen atom, there are quadrupo- 
lar contributions coming from the currents flowing 
along the neighboring elementary cells which are 
not pointing out to this atom (the intensities of the 
currents are different in the two neighboring cells) ; 
see Fig 15] b) numerically, we get a value of 50 [Oe] 
for the field. 

• for all other atoms (except the on-rung oxygen) 
we get octupolar contributions (as shown in FigH 
a)) coming from currents flowing through unit cells 
other than that of the atom: one piece stems from 
the single octupole on the leg opposite to where site 



x resides and the other from a sum of double oc- 
tupoles along the ladder (a series oscillating with a 
(2kFo)^ 1 periodicity and decreasing as 1/r 5 ); nu- 
merically, we get a value of 10-=-20[Oe] for the field. 

We note that both types of contributions produce a stag- 
gered field pointing in a direction perpendicular to the 
plane of the ladder, which possesses the periodicity of 
the o-OCP. The dipolar component and one of the terms 
in the octupolar series are illustrated in the figure below. 



O O 





b) 



o c> o 




FIG. 2: The two main contributions to the amplitude of the 
magnetic field at lattice site x (orange circle) as a result of 
the o-OCP flow: (a) the octupole moments (single octupole 
= green square) which should be summed up for all elemen- 
tary unit (squares) along the ladder (except the one which 
contains lattice site x) and (b) the quadrupole magnetic mo- 
ment that is present only if x is one of the on-leg oxygens; 
the field comes from the asymmetry of OCP amplitudes in 
two neighboring unit cells (non-zero contribution indicated 
in red); crosses (dots) correspond to single magnetic dipoles 
pointing into (out of) the page 

This local magnetic field B(x) break SU(2) spin- 
rotational symmetry. Usually, such symmetry breaking 
may occur in the presence of spin-orbit, which adds a 
L e ff(x) ■ S(x) term to our Hamiltonian (it should be 
emphasized that both vectors are taken at the same 
point). Tools required to deal with such a pertur- 
bation are readily available. We simply need to tai- 
lor these to our case. Both operators are two o-band 
fcrmion combinations, which, when coupled, will give 
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four fermion products (to first order). Detailed analy- 
sis shows that such L e tt{x) ■ S(x) terms indeed cover 
the presence of magnetic field provided that hypothet- 
ical (periodic along x-axis) held B(x) would have had 
only z-component. For example we may check the form 
(in bosonization language) of four fermion operators de- 
scribing the coupling of the SDW with the spin degree of 
freedom of carriers (an example of this type of coupling 
is: ipl^ix^oB^^i^lL/R^i^^oL/R^ix))) Thus, one 
can use the theory of spin orbit coupling (as a pertur- 
bation) to solve our problem. The z-axis orientation of 
our perturbation is a consequence of the high symmetry 
of the lattice, a point which we emphasize in the next 
section. 



C. Side effect of the o-OCP: new terms in the 
hamiltonian 

1. First order term: spin dependent hopping 

Because of the local magnetic fields that are generated 
by the o-OCP, the on-site energies of each atom (except 
for the on-rung oxygens), will be different for carriers 
with spin up and spin down. This can be viewed as local 
Zeeman splittings, which also affect overlap integrals for 
different spin rotations. The d-p orbital overlaps, and 
thus the hopping amplitudes, will depend on S z . This 
suggests that spin-dependent hopping will be present. 
The additional hopping term is defined by t acr i = it ■ a, 
where a are Pauli matrices. Our argument implies that t 
has only z-components. In general there can be also other 
sources of t like coupling of electrons orbital momentum 
L with magnetic moment of each triangle. 

For Cu-O lattices, we can make generic statements re- 
garding the possible spin-dependent hoppings caused by 
interactions with finite out of plane moments (the treat- 
ment is similar to a crystal field splitting with a specific 
direction in space): 



(<1 



x*—y- 



\L\d n ) 



(12) 



Here d x i_ y i is the ground state d-orbital of the Cu atom 
located at some site i; \d n ) is an excited d state of that 
atom; t p d n is the hopping amplitude between the Cu 
atom in the state \d n ) and a neighboring O atom in 

one of the p orbitals represented by the state |p 7 ); L 
is the orbital moment. One has to remember that sym- 
metry requires some of the hopping integrals to vanish 
(for instance t Pj d xz = and similarly L\d^ z 2_ r %) = 0). 
Thus lattice symmetry implies that t has always only a 
z-component. 

Taking this extra term into consideration, first order 
effects can be easily incorporated in our tight binding 
model: the problem can still be factorized in the spin 
variables (t has only a z component). The magnetic field 



has o-band wave-function symmetry, thus Zeeman split- 
tings (and also overlap integrals variations) will have the 
same intra-cell symmetry. For the OCP case, the split- 
ting is purely vertical and it oscillates with kp Q period- 
icity from unit cell to unit cell (whereas in the more fre- 
quently discussed Rashba case it is horizontal and the 
same for all the unit cells), which gives constructive inter- 
ference within that band. The result is that the o-band 
is split - the hopping parameters of spin up and down 
are slightly different. 

There are three consequences of this type of splitting: 

• the Fermi velocities of spin up and down states are 
not equal anymore 

• there is a mismatch of their Fermi vectors 

• the two bands may contain different amounts of d 
and p orbitals 

We do not need to discuss the last item, because the 
eigenvectors of bands pairs are very similar. Let us in- 
stead focus on the first two. The difference in the veloci- 
ties, translated into boson field language, yields terms of 
the form c +II s +, which affect the values of the param- 
eters K c ^ s+ and the form of correlation functions. The 
modification of K CtS+ has to be taken into account in the 
initial condition of the RG flow; the impact on correla- 
tion functions will be described in detail in Sec. IIV1 The 
Fermi vector mismatch has an adverse impact on perpen- 
dicular (OCP case) spin scattering, because momentum 
is not conserved during that process, any longer. This 
effect will be accounted for much in the same way as 
one treats slight departures from commensurability with 
doping, in Mott transitions. 

In our case, we may estimate \i\/t < 10~ 4 (considering 
that the magnetic moment of a single loop is A ~ O.lfig). 
This value gives 5kFi/kFi of the same order, so all these 
effects can indeed be treated as a perturbations. 



2. Second order terms: spin-flip scattering channels 

Let us turn to the next order in perturbation, which 
yields four fermion operators. In general these can be 
written in a Dzyaloshinskii-Moryia (D-M) form 



D-M 



(13) 



The value of the D-M interaction between two Cu sites 
i and j separated by an O site I can be computed from the 
third order spin dependent hopping perturbation. The 
general expression is (where in the limit of small interac- 
tions e is equal to the difference in the O and Cu on-site 
energies) 

Dij = — ^(Uitij+tuti^^ituUj-tjitu+tiixtij) (14) 



The non-zero value of comes from spatial inho- 
mogeneity of (by its definition periodic) OCP pattern, 
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which implies ty 7^ tu. Taking into account t ~ 10~ 6 eV 
and using the fact that e < t in Cw-based ladders (and 
when U —> as for the C2S2 phase, e determines inter- 
level distances), we can estimate D ~ 10~ 5 eU (it is of 
order ~ t(j) 3 ). This value is one order of magnitude 
larger than magnetic dipole-dipole interactions expected 
for such compounds (see Appendix). 

We see that in the OCP case, where t\\z, the cross-term 
term in the second parenthesis of Eq. (|14[) is zero, and 
hence D is perpendicular to the plane of the ladder. It 
also implies that this four-fermion operator will be pro- 
portional to the asymmetry between neighboring Cu-Cu 
links (il-th and Ij-th). We note that if the local magnetic 
moments were constant along the ladder, D = 0. 

Thus, the periodicity of the symmetry breaking per- 
turbation plays a critical role here. The OCP is super- 
imposed on a DW with a (2kp )^ 1 periodicity in real 
space, and the Fourier transform of D(x) - which origi- 
nates from the local magnetic fields - has components at 
\q\ = 2kp only, i.e corresponds to a backward scatter- 
ing vertex. This mirrors the fact that the spin-rotational 
symmetry breaking physics brought about by currents is 
due to the o-band and can only involve carriers scattering 
in this band. 

When Dij\\z one may deduce that the processes which 
explicitly break spin-rotational symmetry are of the form 
S x S y . Using spin- flip operators we have terms like 
S^S- + h.c Physically this indicates that spin is not 
a conserved quantity anymore and that new instabilities 
may arise as a result of the additional scattering chan- 
nels. The new term in the hamiltonian is 



H D -M = 9fo / IpLoa^Roa^LoalpRoa 



(15) 



where gf is the 2kp component of the Fourier transform 
of Dij . This can be bosonized in the standard way: 



H D -M = 9fo J drcos(28 so ) = g fo J dr<. 



Since spin-flip terms are involved in the fermionic opera- 
tors, none of these will give simple density-density contri- 
butions to the bosonized expression of the hamiltonian; 
instead they produce non-linear interaction terms which 
require RG perturbative treatment. They may thus open 
a gap in the spectrum and cause new types of orderings 
in the ladder (which were previously forbidden as they 
explicitly break SU(2) symmetry). We investigate this 
possibility in the next section. 



III. ANALYSIS OF THE SU(2) SYMMETRY 
BREAKING TERMS USING BOSONIZATION 
AND RG ARGUMENTS 

The standard way to deal with four fermion operators, 
which produce cosine terms in bosonization language, is 
to follow an RG approach. This was extensively discussed 



in the context of SU(2) symmetric two-leg ladders (see 
for example Sec. Ill in |29j) and we transpose it to the 
present situation. 

As discussed in the previous section, we add the gf a ~ 
D(2kp )-type contributions to the Hubbard-type inter- 
actions and, following the standard procedure for cosine 
terms in the bosonization representation, we perform an 
RG analysis of the set of four fermions operators. Be- 
cause spin-flip processes are now allowed, we have to add 
scattering events where - for instance - two incoming 
fermions are in a spin-up state and the two outgoing par- 
ticles are in the spin-down state. 



A. The extended system of RG equations 

We proceed along the same lines as in our previous 
studies of two-leg Cu-O ladders; however additional sub- 
tle factors have to be taken into consideration. The fact 
that D(q) has always just one component in momentum 
space simplifies the analysis and enables us to exclude 
some processes 

One should take the system of RG equations derived 
for two-leg Cu-O Hubbard ladders ( Eq. [17] to Eq. [31] 
in Ref. |291 ) and add to it two equations describing the 
renormalization of the gf a terms (intra o-band spin-flip 
scattering) : 



dl 



g fl -(2-(K- 1 +K- 1 ))+ 



P 1 Q 1 (K^ 1 - K^)g f2 (17) 



-g f2 -(2-(K^+K^)) + 
PMK- 1 



dl 



where gf\ 



.5/2 



K^)g n (18) 
and by definition, 



in our case, g^ = initially. We are using the same 
notation as in Ref. 29: K\^ 2 are the LL parameters of 
the eigenmodes in the spin sector, i.e linear combinations 
of K s± such that K\ = P?K S ^ + Q\K S+ and K 2 = 
PiK s+ — Q\K S ^. The initial values of K s ±, P, and Q 
are determined from the Fermi velocities in the two bands 
and interactions of the forward type which are included 
in the LL matrix K (note that Pf + Q\ = 1). Since g/ 
has a cosine form in bosonization language, it does not 
change the initial values of the entries of the LL matrix 

k. 

Basically the new couplings <?/i and gf 2 can be treated 
in a way similar to gi, g 2 ; as they pertain to the 9 +s field, 
these cannot couple with other interactions <?;. The sec- 
ond order g^ Q terms will enter the renormalization equa- 
tions of the LL parameters Ki and cot 2a (eigenbasis an- 
gle) in the spin sector: 



A 



~df 



(9% 



9%) + f(Pi)(9fi9f2 



(19) 
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A ^T = 2 ^ + ~ f ^3fi9f2) (20) 



where, as before: 

/(P 1 ) = (p 1 Q 1 + I^-^i)-i (22) 

/i(Pi) = ((PiQi) 2 + 0.25(P 2 - Q 2 ))- 1 (23) 

We notice that, as one might have expected, only the 
spin sector is affected by the new non-linear terms. The 
modifications that these terms cause to the RG flow (in 
particular as far as the eigenbasis rotation is concerned) 
are not easy to evaluate, because - as we argue below 
- several effects play a role. A qualitative analysis of 
the RG equations only allows one to predict that the 
evolution of cot 2a during the flow will be much slower 
than in the standard case. 

The initial values of the Ki parameters will be different 
from what they are in the standard case, if only because 
of spin-charge mixing (following our discussion in Sec. 
Ill C f 1 spin up and down are not degenerate any longer). 
This effect was partially accounted for by Moroz et al<^: 
they found that it can be implemented through a proper 
change of the scaling dimension of cos(</>) (Eq. [38] in 134 ). 
More explicitly, K s+ and K c+ will mix (just as K s ± do, 
in response to particle-hole symmetry breaking) with a 
mixing parameter proportional to e = (Vp-j- — Vf\)/Vf 

Assuming that the o-band is decoupled from the 7r- 
band, we can easily carry over Moroz's results to the 
ladder case and conclude that the mixing of spin and 
charge generated by the spin-orbit coupling is irrelevant 

v° 

in the o-OCP case, because here e <C l crit — p^- 

Yet, another effect does play a role. As we pointed out 
in the band structure discussion, because of the Fermi 
vector mismatch, perpendicular spin scattering is re- 
duced compared to parallel spin scattering (as is the case 
with umklapp scattering, slightly away from fci? = 7r/2). 
This causes a shift of the phase boundaries in the ground 
state phase diagram, but it does not generate any new 
phases. To demonstrate this explicitly requires a full 
blown RG treatment but heuristic arguments pertaining 
to the RG flow in the C2S2 phase - presented in the next 
section - give clear indications of this trend. 

B. Relevance of the gf terms 

The overall structure of the phase diagram that 
emerges from the RG analysis, when these additional 
terms are included, is similar to that obtained in our pre- 
vious papers, so we discuss here the stability of the phases 



that were found for ladders possessing SU(2) symmetry. 
Let us start from the half filled case when there is on av- 
erage one carrier per copper atom and move towards the 
highly doped regime when the 7r-band is nearly empty. 

The low doping C1S0 phase is not affected, since the 
o-OCP does not exist, and hence D{2kp ) = 0. How- 
ever, to be on the safe side, we have checked that even 
if D(2kp ) 7^ interactions of the gf Q type are irrele- 
vant. This can be justified by symmetry considerations 
regarding the eigenbasis: in the low doping regime 
is the proper basis of LL modes, thus a small g/ , which 
induce instability in the B 07T basis cannot be relevant in 
this doping range. It must be interband spin flip scat- 
tering which causes a relevant perturbation of the CI SO 
phase - the possibility of the other types of spin flip 
interactions are briefly discussed in the Appendix. 

A significant difference, caused by gj Ql is found for the 
massless C2S2 phase that was present in a finite, inter- 
mediate doping range of the SU(2) symmetric problem. 
We recall, that for this phase, P ow is the proper basis 
of the spin sector, thus interband instabilities are sup- 
pressed. Adding gf a favors a flow towards this basis. We 
can simplify our analysis even further if we assume that 
in given range of dopings, the RG equations describing 
the evolution of the o-band are decoupled from the rest 
of RG system. This allows us to work with Eq. (| 1 7[) and 
the equation (derived previously 2 ^) describing intraband 
backscattering: 

^=g-(2-{K 2 + K,)) + P\Qi{K 2 - Kjg 

g a was proved to be irrelevant since the initial K 2 > 1 
and we would need a large value of Q\ in order to open 
a gap in <j) so . Roughly speaking K so was approaching its 
fixed point K so — > 1 too slowly and the flow was finally 
only marginal. In our extended case gj tends to increase 
K so or at least (in the beginning) protects its constant 
value. The arguments given in 29 work exactly in the 
opposite way for the cos 9 so instability, which is relevant 
provided K so > 1. This can be quantified by studying 
an invariant of our extended RG system: 

A 2 = (K S0 -l) 2 -(g 2 + g%) (24) 

We see that by adding the new terms A 2 is shifted (from 
A a; 0) to negative values; it is known that an imagi- 
nary A in the Kosterlitz-Thouless RG flow corresponds 
to a divergent trajectory that is to a relevant instability 
destroying the LL fixed point. The only possibility cor- 
responding to that case is that gf becomes the relevant 
perturbation. 

Our earlier observation that momentum mismatch re- 
duces the initial value of the perpendicular spin scatter- 
ing interaction while leaving those of the parallel and 
spin flip processes unchanged reinforces the inequality 
(K so — 1) > g±. We then study a Kosterlitz-Thouless 
RG flow with all the g± terms smaller than the g» a terms; 
in that case we find that we move away from the SU(2) 



separatrix and scale towards larger K SOl which favors the 
new intermediate phase (with the 9 field ordered). 

One finds that a new gap opens in the spectrum of the 
9 so field. We find that C2S2 -> C2SV, with an order- 
ing within the o-band. Using a bosonic version of order 
operators we find that SDW in the plane of the ladder 
(SDW X in the o-band) is favored. The order operator is 
defined as follows: 

O x sdwo(x) = ^Ir^oLi(x) + ^l^oLtix) 

~ exp(— 2ikF x) exp(zv / 2</'co(^)) cos(V29 so (x)) 

The interaction breaking spin-rotational SU(2) symmetry 
dominates in this regime, so it is clear that the resulting 
phase also breaks this symmetry. 

This simplified analysis is confirmed when we perform 
the actual RG calculation including the additional terms. 
Furthermore, we are able to obtain the exact boundary 
between the gf D and g dominated phases (i.e the bound- 
ary between the SU(2) symmetric and symmetry broken 
phases). We also see that the assumption K so ~ const 
is justified, so we can approximate the value of a gap by 
taking the following RG flow: ^ = 0, from which one 
deduces that a proper estimate of the gap value is 3 -!: 



A o5/Y 



(25) 



Assuming K so < 2 (strong, but local, Hubbard inter- 
action limit) we find A ~ 0.01D, which gives quite a 
small but experimentally accessible transition tempera- 
ture G (1, 10)mK. 

For the C2S1 phase the additional gf terms compete 
with the interaction g a = gi + <?2, so the theory of spin flip 
effects in the case of a single chain applies. The reason- 
ing is very similar to that for the C2S2 phase, except 
that now g dominates. To put it on a more quantita- 
tive level: as long as the initial gj is smaller than g a — 
which is a reasonable assumption, because g a is a rele- 
vant perturbation dominant in this doping range - the 
ordering of the (f> so field overwhelms spin-flip processes. 
Giamarchi and Schulz have found that the critical line 
for this problem is at g/ = g Q ^ . This assumption can 
fail only close to the low doping boundary of the C2S1 
phase, since the strength of g decreases as one reaches 
5 C 2- We can thus conclude that, except for a shift of the 
lower boundary, the CDW state is unaffected by the spin- 
orbit coupling. The new phase diagram, in the presence 
of gf is summarized in the figure below (Figj3]). 



IV. DISCUSSION 

The breaking of SU(2) symmetry, that occurs in the 
presence of the OCP, opens a gap in the spin mode of 
the o-band. In the range of dopings where previously a 
critical line was found, the C2S2 phase is replaced by 
a different ground state, namely a phase with in-plane 
SDW X denoted by C2S1'. An important question is 
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FIG. 3: Modified phase diagram (to be compared with Figfj]) 
of the two-leg Cu-O ladder when intra o-band spin-flip pro- 
cesses are included: as before S is the hole doping and t pp 
inter-oxygen hopping. In the intermediate doping range - 
where a C2S2 phase was previously found in the SU(2) sym- 
metric problem, a C2S1 ' state is present, which has a SDW 
in the plane of the ladder as a result of the opening of a spin 
gap. 



whether this unusual phase is detectable, or if any small 
symmetry breaking will make it undistinguishable from 
the standard, high doping C2S1 phase. To answer this 
question we present a discussion along three lines: the 
first one is devoted to high temperature correlation func- 
tions, the second to transport properties and the third to 
magnetic properties. A limitation of the discussion pre- 
sented below, which affects experiments, is the energy 
scale: the effects will be clearly visible only in a temper- 
atures range below the gap A. 

One may also wonder whether the existence of the 
OCP is the only source of SU(2) symmetry breaking. As 
the Cu-O lattice forms a fairly rigid structure, it is stable 
both against crystal field d levels splitting and standard 
L — S coupling ~ Z 2 on heavy Cu ions (with Z = 29), 
and, insofar as we may neglect size effects, the presence 
of an OCP is the source of symmetry breaking. A de- 
tailed discussion of these issues will be presented in the 
Appendix. 



A. Correlation functions for the LL phase 

Correlation functions for LL with spin-orbit interac- 
tions have been evaluated both for two- 3 "* and four- 
fermion 36 operators. Let us briefly summarize here their 
main features. Chiral separation is still preserved and 
correlators can still be factorized in terms of two modes, 
but not spin and charge anymore. Interactions and spin 
orbit coupling (L-S) reinforce one another; in the pres- 
ence of L-S, the density of states decays faster than for 
free fermions, and the difference between the velocity of 
the two modes SV = V\ — is proportional to t. Ac- 
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cording to2£ the exponents are slightly changed when a 
term proportional to a difference of the velocities is added 
(we use these results slightly modified, because it was 
Rashba-type spin-orbit coupling treated in Ref. [3(| while 
we are interested in Zeeman-type coupling); for example 
the behavior of our dominant fluctuation is given by (at 
zero temperature, up to first order in 8V): 

O x SDWo (r) ~ cos(2fc i , o2 ;)r-^+ 1 /^°- r ) (26) 

where 6 s ~ SV. As a consequence the four fermion op- 
erators have a non-zero conformal spin proportional to 
SV. 

In our case intra- o-band L-S coupling is quite weak 
(this small value enabled us to neglect some effects of 
spin-charge mixing during the RG flow). Accordingly, in 
the high temperature LL phase, the functional (power 
law) form of the correlators will be the same for the 
C2S1 and C2S1 ' phases, and the only difference comes 
from small changes in the exponents. Taking into ac- 
count small (and oscillating) value of 5k (SV), and ex- 
perimental difficulty with extracting exact value of bare 
LL exponent Ki Q , we doubt that C2S1 ' existence can be 
shown in this way. 

More visible signatures can be expected at low temper- 
ature where the 9 so field is ordered, as all the observables 
depending on magnetic susceptibilities - such as NMR 
responses - are affected (notwithstanding the relatively 
smallness of the gap) : if the 6 so field is locked, then av- 
erages over the (f> so field go abruptly to zero. According 
to2£ one expects quite abrupt suppression of the o-band 
electronic part of the Knight shift K and relaxation rate 
(1/Ti) at T < A (for example K ~ (cos^ so )). 

B. Transport properties: analysis of the impurity 
problem 

Let us briefly discuss here what happens if a non- 
magnetic impurity is introduced in the two- leg Cu-0 
ladder. This problem was discussed in detail recently^; 
here we focus on the consequences of spin-rotational sym- 
metry breaking. Specifically, we investigate the doping 
range where the C2S2 phase was stable in the absence 
of spin-orbit scattering. In the previous section we have 
shown that K so behaves differently at the lowest tem- 
peratures, because of a gap opening due to spin-orbit 
terms. So the temperature dependence of the transmis- 
sion through a single impurity - which is connected to 
the renormalized backscattering V [l;K] - will also be 
affected. 

Spin-orbit coupling is present only inside the o-band, 
so we may restrict our analysis to reflection (transmis- 
sion) of carriers inside this band. In such a simplified 
case, the RG equation for backscattering on the impu- 
rity is given by 

^- = {2-K so - K C0 )V (27) 



instead of Eq. [15a] in Ref. 37, where it was assumed 
that K so = 1). From this, it is clear that increasing K so 
(which promotes an ordering of the SO field) can cause 
V to become irrelevant. More precisely, this operator 
is always relevant at the beginning of its flow, because 
K co < 1 (and this value is constant in the absence of 
umklapp terms), but this trend may change as the en- 
ergy is decreased. This would give rise to an unusual 
temperature dependence of the conductivity. 

On the other hand, if the initial value of K co is suffi- 
ciently small (strongly interacting case), we would arrive 
at an open boundary problem (V relevant). In our pre- 
vious paper—, intraband Kondo physics was suggested 
for this case. The following RG equation describes this 
regime (assuming open boundary conditions (OBC)) 

^ = J (l K^) (28) 

We see here that if SU(2) symmetry is broken by the o- 
OCP, this coupling becomes relevant, which confirms our 
previous predictions. In the limit of quite strong Hub- 
bard interactions K co — > 1/2 - which is the case when 
K so 3> 1 - it can give rise to the unusual regime predicted 
many years ago by Furusaki and Nagaosa^ who took 
second order tunneling processes into account, showing 
that double spinon tunneling events may dominate the 
physics around the impurity. The competition between 
OBC Kondo physics and spinon transmission will take 
place in the strong interaction limit (the barrier asym- 
metry will also play a significant role) , but this is beyond 
the scope of the present paper. 

An experimental consequence of this effect, is that in 
the phase where spin-orbit coupling dominates, the usual 
line broadening caused by impurities that is expected in 
an NMR measurement, will not be seen. This is because 
spinons can be transmitted through non-magnetic impu- 
rities (and hence the backscattering in the spin sector is 
expected to be weak). 

C. Magnetic properties 

Obviously, the local magnetic field will affect the 
Knight shifts of the NMR lines. The magnetic field is 
alternating and incommensurate with the lattice: as a 
result of such local magnetic field variation NMR lines 
will not have standard (thermal broadened) shape, but 
it will depend on the number of atoms feeling a given 
value of the field N(Bi oc ). The incommensurability with 
the lattice allows us to rephrase this question and ask 
what is the density of states for a cosine curve. Thus the 
NMR line-shape resemble ID density of states with two 
peaks line-shape (ID density of states has 1/tj singularity 
coming from the bottom of the cosine band) . The split- 
ting is equal to 2B max and, because the effective field is 
different at each atom inside unit cell, thus different line- 
shapes are expected: the widest on on-leg oxygen atom, 
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and non-perturbed (zero split) on the on-rung, central 
oxygen. This effect should be distinguished from the one 
described in the previous paragraph, in total OCP mod- 
ifies NMR lines in two ways: it suppress the amplitude 
of the satellites, of the impurity at the origin, and splits 
the central peak. Unfortunately estimated values of this 
splitting are at the limit of the current experimental res- 
olution, so this primary effect might be difficult to de- 
tect. One can alternatively detect anisotropic g-factors 
through Zeeman effect. 

There will also be a signal coming from the magnetic 
ordering itself. We predict that there is a magnetic mo- 
ment component perpendicular to the ladder plane due to 
B e ff induced by the o-OCP (primary effect), but there is 
also a secondary effect: from our RG procedure we have 
found an in-plane component corresponding to a SDW X , 
with the same periodicity as for the staggered moments 
along z. The fact that both components have the same 
periodicity in k space, 2kp 0l is quite important: they 
can add up creating a moment tilted out of plane. The 
tilting angle will be determined by their relative ampli- 
tudes. We have started with \x z = 0.1/is which induce 
an ordering gap A. The value of /j, x is connected with 
the average \x x ~ (cos SO ) so it will be reduced with 
increasing temperature. Even at zero temperature quan- 
tum fluctuations cause (cosf? so ) 7^ 1. In order to estimate 
the average at zero temperature we compare the values of 
the gap we found in Sec IIIII with the one found when the 
LL velocity (kinetic energy) goes to zero Aqo . It is known 
that Aqo ~ \D\. The ratio A/Aoo = 0.01 is proportional 
to (cos# so ) and allows us to estimate ji x w 0.01/x_b, from 
which we deduce a tilting angle of the order of •& ~ 6deg 
(counted from vertical axis). This kind of modulated 
(in space) magnetic moments composition should be de- 
tectable using elastic neutron scattering, but quite low 
temperature (~ lOmK) would be required for this ex- 
periment. 



V. CONCLUSIONS 

The main question which we addressed in this work 
was: what is the feedback effect of an OCP on the prop- 
erties of the ladder, if one takes into account the coupling 
between the electrons and the induced orbital moment? 
Such feedback can in fact potentially provide an addi- 
tional method to probe the existence of the OCP; it does 
not require a direct measurement of the induced - and 
quite tiny- magnetic field. Detecting such feedback can 
however prove challenging, due to two important con- 
straints. First, one would need to obtain highly doped 
(~ 25%) two-leg ladders; second, experiment have to be 
performed at temperatures below A. Still, we hope that 
our present detailed insight into the low energy proper- 
ties of the intermediate phase will invite experimental 
and numerical investigations. 

A non-trivial issue is how reliable our value of A is. We 
have used a weak coupling RG approach which is reason- 



able in the case of the otherwise gapless C2S2 phase: the 
statement A < \D\ is justified. Thus, in our framework 
the only way to increase A, assuming a given amplitude 
of the OCP (~ fi z ) is by taking larger A coefficient in 
Eq. (fT2|) (this would give larger |tj). This may arise be- 
cause of particular, unknown property of copper oxides 
and/or constructive interference of L-S couplings defined 
on the top of the DW. Larger gap value will push the tilt 
angle d towards larger values. 

Beyond the experimental manifestations, which we 
have discussed in Sec lIVl the magnetic moment pattern 
should be detectable using elastic neutron scattering. 
Such experiments were recently performed in the pseu- 
dogap phase of 2D HTS C 12 i 15 . The data was interpreted 
as evidence for out-of plane magnetic moments with an 
alternating tilt. There is a similarity with our prediction, 
but one must not forget that there are several important 
differences between the OCP that we have discussed and 
the one that could serve to interpret the experiments in 
cuprates: in our ID case the pattern is modulated in 
space with a kFo periodicity (this modulation was in fact 
critical to obtain a non-zero coupling with the spin de- 
gree of freedom) and has a different internal symmetry: 
6*i instead of 62 in Varma's notation^. One must keep 
in mind these caveats, nonetheless our result invites one 
to look for similar effects in two-leg ladders. 

To summarize: despite certain similarities, the 2D 
cuprates studied in Ref. [T2IIT5I cannot be treated as an 
experimental realization of our problem. A good candi- 
date would be a doped ladder, such as the so called tele- 
phone compound {Srn- x Ca x Cu2iOi\) 1 provided suffi- 
cient doping levels can be attained (from the standpoint 
of solid-state chemistry, this is a highly non-trivial task) . 
Neutron diffraction or /iSR experiments would be desir- 
able for such a ladder system. 

An important point is the possible role of apical oxy- 
gens. In the framework of a tight binding model the 
new hoppings will make the two bands (o and 7r) more 
asymmetric; this implies larger effective t pp , which shifts 
the C2S2 phase towards smaller dopings 6 (see FigfT]). 
New conducting paths will also increase the OCP am- 
plitude for a given value of t pp , shifting the critical t™ n 
to smaller values. As a result the new phase C2S1 ' will 
be easier to access experimentally in the (5 — ^2.) plane. 
The new current paths can also give rise to higher order 
L-S couplings: for example pi can induce additional cur- 
rents flowing in the loops involving the apical oxygen. We 
may conclude that including apical oxygen should favor 
the effects predicted in this paper. 

This statement, about a positive role of apical oxygens, 
is in qualitative agreement with recent work^i studying 
currents emerging because of particular spin-textures. 

A non-trivial outcome of this work is that it under- 
scores the hidden connection between a phase with OCP 
and a pseudo (partial) gap in the spin sector. The emer- 
gence of such a gap shields the intermediate, previously 
marginal, phase from weak perturbations such as inter- 
ladder hopping or disorder. In conclusion our work shows 
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that a strong perturbation is caused by the presence of 
the OCP and validates the phase diagram for two leg 
Hubbard ladders which had been obtained in Ref. [30|. 

VI. APPENDIX: OTHER SU(2) SYMMETRY 
BREAKING MECHANISMS 

In this Appendix, for the sake of completeness, we 
briefly discuss alternative sources of spin-rotational sym- 
metry breaking. A few of them come readily to mind: 
dipole-dipole interactions, standard L ■ S spin-orbit cou- 
pling and Rashba mechanism. 

1) Direct spin-dipolar interactions. The value of this 
spin anisotropy depends on the average distance between 
carriers as ~ l/(r) 3 . The (r) is obviously changing when 
one add carriers into the ladder which implies a doping 
dependence. We can use values found for the half filled 
case in a similar compound^ and simply assume that 
carriers are more diluted. For the most interesting phase 
(C2S2) there are at least two times less carriers which 
gives reduction by a factor 8, thus we have estimated the 
dipole-dipole anisotropy to be one order of magnitude 
smaller than OCP effects. 

Besides, in the C2S2 regime, OCP effects are predom- 
inant because they possess the proper o-band symmetry, 
in contrast with the generic form of dipolar interactions. 
Thus we assume that this mechanism cannot compete 
nor mimic the OCP effects predicted in this paper. 

2) The spin orbit L- S coupling is a fairly typical effect 
for heavy atoms. Thus we should check whether it could 
induce much stronger effects than those predicted in our 
analysis. 

The key point is the symmetry of the lattice: so long as 
copper atoms are sitting in a C4 symmetric environment, 
hoppings to and from all the neighboring oxygens will 
conspire to give a zero net outcome. We can also invert 
this argument and say that if electron ordering finds a 
way to lower the rotational symmetry of a lattice then it 
would provide an efficient way to increase the A coefficient 
in Eg IT21 and thus enhance effects discussed in this paper. 
This property does not hold if the lattice is distorted (in 
the Peierls or Jahn- Teller cases for instance) or warped. 
As the lattice is quite rigid (and we are not working at 
half filling) for ladders, these effects are not expected to 
be significant in the bulk. The situation may be different 
at an interface, but we are not dealing with such a case in 
this paper. This is the reason why we focused on OCP- 
induced effects. 

The only exception might be if the primary OCP in- 
volved apical oxygens. Let us make a gedanken experi- 
ment and assume that there is a current flowing through 
some of bonds toward these atoms. This might enhance 
the total amplitude of the OCP, lowering the value ipp"*' . 
It can also be a higher order perturbation induced by the 
in-plane SDW(x). The conducting bonds tend to be 
shorter which can displace apical oxygens. The emer- 
gence of such distortion will lower the lattice symmetry 



from C4 to C2 allowing for the non-zero net expecta- 
tion value of the spin-orbit coupling on Cu atoms. This 
would cause a S z dependent hopping, as can be seen from 
Eq. |12| which shows that the displacement is equivalent 
to admixing d x 2_ y 2 orbitals with the apical oxygen band. 
Importantly, this effect is caused by the OCP, thus it 
should have the same periodicity and phase. This ad- 
ditional L-S coupling enhances the OCP induced effects 
discussed above. 

2) Another possible mechanism of interest in the 
cuprate case, giving spin-dependent hopping, arises from 
the Rashba effect which couples directly electron motions 
to their spins. It describes the precession of carriers in 
the presence of a strong electric field. The Rashba Hamil- 
tonian is: 

H R =a(kxn)-a (29) 

where a is a coupling constant and ft is a unit vector 
in the electric field direction. Obviously, this coupling 
(which depends on the total momentum of the carrier) 
will be different for the o and tt bands, because of the 
inequality kp a > fcpv (one should not forget that /cpv 
has large on-rung components). 

Firstly, let us consider the bulk of a ladder material. 
The electric field which we expect will lie in the plane of 
the ladders in a direction perpendicular to the legs. It 
will result from the arrangements of DW in neighboring 
ladders (so again it will be restricted to the high dop- 
ing regime, and again have a {2kp Q )~ l periodicity in real 
space). There are two limitations here: these arrange- 
ments can only be stable for certain dopings (quarter- 
filling) and they are expected to be quite small in ampli- 
tude. 

Straightforward analysis (with again crucial element 
of DW modulation) of Eq. [55] gives a spin-dependent 
hopping parameter t ~ [0, 0,t'], so again it has only a z- 
component. It implies that the analysis done in the core 
of this paper applies also in this case. Thus, this kind of 
Rashba (if such kind of perpendicular arrangement does 
exist) would enhance the effects predicted before. The 
only difference which would allow us to distinguish this 
case from that induced by an OCP is the magnetization 
component perpendicular to the ladder plane. 

Let us finish with a brief discussion of an interface 
where the Rashba mechanism can induce quite new 
physics. In this case the electric field will be perpen- 
dicular to the ladder plane, because it arises from charge 
polarization effects at the interface. On the interface, 
when charge differences are very high (polarization catas- 
trophe) at very small (interatomic) distances, this effect 
can be really huge. What is more, in contradistinction 
with the mechanisms that we discussed until now, the 
amplitude is not alternating in space (the sign of the 
coupling is not oscillating). 

In the case of strong electric field there are only non- 
diagonal (spin-flip) hoppings; level splitting is induced by 
the Rashba term Eq. [5pJ (if the electric field E is in the 
z direction, Hr can be expressed as ~ [t', 0, 0] • a in the 
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o-band and approximately by [0, t' , 0] • a in the 7r-band). 
In this case spin and charge are strongly mixed thus spin- 
up and spin-down are not a good quantum number any 
longer. 

Also spin-flip interactions will be highly non-trivial in 
such a model: we may find t%j with perpendicular ori- 
entations for instance when the two scattering carriers 
are from two different bands. Thus it is the last term in 



the parenthesis of Eq. (|14p which is non-zero; we are then 
dealing with interband scattering. We note that both t he 
in the ladder plane, so that, once again, -\\z (spin-flip 
process). We expect very different physics emerging in 
this case, which is out of scope of this paper. A detailed 
RG analysis will be postponed to a future publication 
devoted to surface effects. 
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